Discontinuous Euler instability in nanoelectromechanical systems 
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We investigate nanoelectromechanical systems near mechanical instabilities. We show that quite 
generally, the interaction between the electronic and the vibronic degrees of freedom can be ac- 
counted for essentially exactly when the instabihty is continuous. We apply our general framework 
to the Euler buckling instability and find that the interaction between electronic and vibronic de- 
grees of freedom qualitatively affects the mechanical instability, turning it into a discontinuous one 
in close analogy with tricritical points in the Landau theory of phase transitions. 
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I. INTRODUTION 

The buckling of an elastic rod by a longitudinal com- 
pression force F applied to its two ends constitutes the 
paradigm of a mechanical instability, called buckling 
instability.Ii' It was first studied by Euler in 1744 while in- 
vestigating the maximal load that a column can sustain.'^' 
As long as F stays below a critical force Fc, the rod re- 
mains straight, while for F > Fc it buckles, as sketched in 
Fig. [T^-b. The transition between the two states is con- 
tinuous and the frequency of the fundamental bending 
mode vanishes at the instability. 

There has been much recent interest in exploring buck- 
ling instabilities in nanomechanical systems. In the 
quest to understand the remarkable mechanical prop- 
erties of nanotubes,'^2lEI there have been observations of 
compressive buckling instabilities in this system.'^' The 
Euler buckling instability has been observed in Si02 
nanobeams and shown to obey continuum elasticity 
theoryF^ There are also close relations with the recently 
observed wrinkling and possibly with the rippling of 
suspended graphene samples. Theoretical works have 
studied the quantum properties of nanobeams near the 
Euler instability,'i2HiIl proposing this system to explore 
zero-point fluctuations of a mechanical mod^^^^or to serve 
as a mechanical qubit.^'^ 

In this work, we study the interaction of current flow 
with the vibrational motion near such continuous me- 
chanical instabilities which constitutes a fundamental 
issue of nanoelectromechanics.'^ Remarkably, we find 
that under quite general conditions, this problem ad- 
mits an essentially exact solution due to the continu- 
ity of the instability and the consequent vanishing of 
the vibronic frequency at the transition ( "critical slow- 
ing down"). In fact, the vanishing of the frequency 
implies that the mechanical motion becomes slow com- 
pared to the electronic dynamics and an appropriate non- 
equilibrium Born-Oppenheimer (NEBO) approximation 
becomes asymptotically exact near the transition. Here, 




FIG. 1: (color online) Sketch of a nanobeam (a) in the flat 
state and (b) the buckled state with two equivalent metastable 
positions of the rod (solid and dashed lines). An equivalent 
circuit of the embedded SET is shown in (c). 



we illustrate our general framework by applying it to the 
nanoelectromechanics of the Euler instability. 

We find that the interplay of electronic transport and 
the mechanical instability causes significant qualitative 
changes both in the nature of the buckling and in the 
transport properties. In leading order, the NEBO ap- 
proximation yields a current-induced conservative force 
acting on the vibronic mode. At this order, our principal 
conclusion is that the coupling to the electronic dynamics 
can change the nature of the buckling instability from a 
continuous to a discontinuous transition which is closely 
analogous to tricritical behavior in the Landau theory 
of phase transitions. Including in addition the fluctua- 
tions of the current-induced force as well as the corre- 
sponding dissipation leads to Langevin dynamics of the 
vibrational mode which becomes important in the vicin- 
ity of the discontinuous transition. Employing the same 
NEBO limit to deduce the electronic current, we find that 
the buckling instability induces a current blockade over 
a wide range of parameters. This is a manifestation of 
the Franck-Condon blockad^^^^ whenever the buckling 
instability remains continuous but is caused by a novel 
tricritical blockade when the instability is discontinuous. 
The emergence of a current blockade in the buckled state 
suggests that our setup could, in principle, serve as a 
mechanically-controlled switching device. 
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II. MODEL 

Close to the Euler instability, the frequency of the fun- 
damental bending mode of the beam approaches zero, 
while all higher modes have a finite frequency.'^' This al- 
lows us to retain only the fundamental mode of amplitude 
X (see Fig. [l^-b) and following previous studies,'^^'^^ we 
reduce the vibrational Hamiltonian tc^®- 



Hvih — 



2m 



-X' 



(1) 



which is closely analogous to the Landau theory of contin- 
uous phase transitions. In Eq. ([T]), P is the momentum 
conjugate to X and m denotes an effective mass. The 
mode frequency uP' ^ \ — F/Fc changes sign when F 
reaches a critical force Fc. Global stability then requires 
a quartic term with a > 0. Thus, for F < F^, {uJ^ > 0), 
X = is the only stable minimum and the beam remains 
straight. For F > F^ (w^ < 0), the beam buckles into 
one of the two minima at ±X-|_ = —vuliP' j a.. 

The vibronic mode of the nanobeam interacts with an 
embedded metallic single-electron transistor (SET), con- 
sisting of a small metallic island coupled to source, drain, 
and gate electrode (Fig. l];). We assume that the SET 
operates in the Coulomb blockade regim^^ and that bias 
and gate voltage are tuned to the vicinity of the conduct- 
ing region between SET states with, say, zero and one 
excess electron. The electronic degrees of freedom cou- 
ple to the vibronic motion through the occupation h of 
excess electrons on the metallic island. Specifically, we 
assume that the electron-vibron coupling does not break 
the underlying parity symmetry of the vibronic dynam- 
ics under X — >■ —X. This follows naturally when the 
coupling emerges from the electron-phonon coupling in- 
trinsic to the nanobeanPni and implies that the coupling 
depends only on even powers of the vibronic mode coor- 
dinate X. The dominant coupling is quadratic in X, 
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X^n, 



(2) 



with a coupling constant g > 0.1^ When there is a signif- 
icant contribution to the electron-vibron coupling origi- 
nating from the electrostatic dot-gate interaction, we en- 
vision a symmetric gate setup consistent with Eq. ([2]). 

In the presence of the vibronic dynamics X{t), the elec- 
tronic occupation n{X, t) of the island is described by the 
Boltzmann-Langevin equatiorPi' 



dn 



= {n, i/vib} + r+(l - n) - r_n + 5J+ - SJ-. (3) 



This equation assumes that the bias is large compared 
to temperature so that tunneling is effectively unidi- 
rectional and the relevant tunneling rates T± for tun- 
neling onto (-I-) and off (— ) the island are given by 
r± = R-\V/2 ± Vg)e{V/2 ± Vg). Here, R denotes 
the tunneling resistances {R ^ h/e^) between island and 
leads, V is the bias voltage, Q{x) denotes the Heaviside 



step function, and we set h = e = 1. Since both gate 
voltage (via the capacitances C and Cg in Fig.[T|:) and vi- 
bronic deformations couple to the excess charge n on the 
island, the effective gate voltage Vg = Vg — gX^/2 com- 
bines the gate potential Vg (measured from the degen- 
eracy point between the states with zero and one excess 
electron) and the vibron-induced shift of the electronic 
energy described by He- The stochastic Poisson nature 
of electronic tunneling is accounted for by including the 
Langevin sources SJ± with correlators (5J-|_(i)(5J+(t')) — 
r+(l - n)S{t ~ t') and {SJ_{t)SJ^{t')) = r^nS{t - t'). 
The vibronic dynamics enters Eq. ([3| through the Pois- 
son bracket {rt,i?vib}- 



III. STABILITY ANALYSIS 

We are now in a position to investigate the influence 
of the electronic dynamics on the vibronic motion. Near 
the instability, the vibrational dynamics becomes slow 
compared to the elec tronic tunneling dynamics. As has 
recently been shown,'22l2n! the effect of the current on 
the vibrational motion can then be described within a 
NEBO approximation in which the vibrational motion is 
subject to a current-induced force —gXn{X, t) originat- 
ing in the electron-vibron interaction ([2]). This current- 
induced force involves both a time-averaged and conser- 
vative force as well as fluctuating and frictional forces, 
resulting in Langevin dynamics of the vibronic degree of 
freedom. 

In lowest order, the Langevin dynamics only involves 
the conservative force which emerges from the average 
occupation ?io(X) in the absence of fluctuations {5J± ~ 
0) and vibronic dynamics ({n, i/vib} — 0). In this limit, 
Eq. (|3| reduces to the usual rate equation of a metallic 
SET so thafi^ 



Vg{x) 



Vg{x) > v/2, 

-f/2 < vg{x) < v/2, 

Vg{x) < ~v/2 



(4) 



with ?; > and Vg{x) 



x^ 12. Here and below. 



we employ dimensionless variables by introducing char- 
acteristic scales Eq — g^ ja of energy, Zq = \fgJot of 
length, and wq = \/ g/m of frequency (or time t) from 
a comparison of the quartic vibron potential in i?vib and 
the electron-vibron coupling He- Specifically, we intro- 
duce the reduced variables x = X/Iq, p = P/mujolo, 



LOot, V = V/Eo, 



Vg/Eo, and r — Rljo/Eq- In 



terms of these variables, we can also write -ffvib + He — 
Eo[p'^/2 + (— e -I- h)x^/2 + x^/4] in terms of a reduced 
compressional force e ~ —muP'/g. 

The current-induced force —xuq^x) has dramatic ef- 
fects on the Euler instability, as follows from a stability 
analysis of the vibrational motion. The (meta) stable po- 
sitions of the nanobeam are obtained by setting the effec- 



tive force fcs{x) 



ex - 



X — xno{x) to zero. Our results 
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FIG. 2: (color online) (Meta)stable (solid blue lines) and unstable (dashed red lines) positions of the nanobeam vs. scaled force 
e for (a) \vg\ < v/2, v < 1/2, (b) \vg\ < v/2, v > 1/2, (c) Vg > v/2, v < 1/2 (for e+ > 1; a similar plot holds for e+ < 1), (d) 
Vg > v/2, V > 1/2, as indicated in the Vg-v plane in (e). The dotted blue line is the result without electron-vibron coupling. 
Notation: e± = 2vg ± f, e = 1/2 + Vg/v, x\ = t, x\ = t — 1, and = (e — e)/(l — l/2t;). Grey indicates conducting regions. 



are summarized in the stability diagrams in Fig. [2j The 
most striking results of this analysis are: (i) The current 
flow renormalizes the critical force required for buckling 
towards larger values, (ii) At low biases, the buckled 
state can appear via a discontinuous transition. 

These results can be understood most directly in terms 
of the potential v^six) associated with /cff(x). Focusing 
on the current-carrying region (shown in grey in Fig. [2] 
and delineated by max{0,e_} < < t+ with e± — 
2vg ± w), we find 

The quadratic term shows that the current indeed stabi- 
lizes the unbuckled state, renormalizing the critical force 
to e = 1/2 + Vg/v when e_ < < e+ (Fig. [2k-b). Re- 
markably, however, the current-induced contribution to 
the quartic term is negative at small and thus desta- 
bilizes the unbuckled state. According to Eq. ([s]), the 
quartic term in the current-induced potential becomes in- 
creasingly significant as the bias voltage v decreases and 
we find that the overall prefactor of the quartic term be- 
comes negative when v < 1/2.*^ It is important to note 
that this does not imply a globally unstable potential 
since the current-induced force contributes only for small 
x^. A sign reversal of the quartic term is also a famil- 
iar occurrence in the Landau theory of tricritical points 
which connect between second- and first-order transition 
lines.'SSl In close analogy, the sign reversal of the quartic 
term in the effective potential ([s]) signals a discontinuous 
Euler instability which reverts to a continuous transition 
at biases u > 1/2 where the prefactor of the quartic term 
remains positive. 

Specifically, when v > 1/2 (Fig. [2|D,d), the current- 
induced potential renormalizes the parameters of the vi- 
bronic Hamiltonian but leaves the quartic term posi- 
tive. This modifies how the position of the minimum 
depends on the applied force in the conducting region 
max{0, £_} < < e+, but the Euler instability remains 
continuous. When v < 1/2, the equilibrium position 
at finite x becomes unstable within the entire current- 
carrying region. This leads to a discontinuous Euler tran- 
sition when e_ < < £4. (Fig. [2k) and to multistability 



in the region < x^ < e+ when e_ > (Fig. [2]:).'2S1 

At the level of the stability analysis, we can also ob- 
tain the current / by evaluating the rate-equation result!^ 
RI{x)/V = 1/4 — [vg{x)/v]'^ at the position of the most 
stable minimum. Corresponding results in the Vg-v 
plane are shown in Fig. [3]a,-f for various values of the 
applied force e. By comparison with the Coulomb di- 
amond in the absence of the electron-vibron coupling 
(dotted lines in Fig.jsj), we see that the Euler instability 
leads to a current blockade over a significant parameter 
range. For v > 1/2, the bl ockade is a manifestation of 
the Franck-Condon blockade,'^^'^ caused by the induced 
linear electron-vibron coupling when expanding Eq. ^ 
about the buckled state. 

In contrast, for w < 1/2, the current blockade is a 
direct consequence of the discontinuous Euler instabil- 




FIG. 3: (color online) Conductance G — RI /V in the Vg-v 
plane for applied force (a) e ^ 0, (b) e — 0.25, (c,g) e — 0.5, 
(d) e = 0.75, (e) e = 1, (f,h) e = 1.25, within (a-f) stability 
analysis and (g-h) full Langevin dynamics {r = — T = 
0.01). Color scale: G = — 5> 1/4 from dark blue to white. 
Dotted lines delineate the Coulomb diamond for g — 0. 
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ity. We have seen above that in this regime, the buckled 
state becomes unstable throughout the entire current- 
carrying region. As a result, the current-induced force 
will always drive the system out of the current-carrying 
region, explaining the current blockade. An intriguing 
feature of this novel tricritical current blockade is the 
curved boundary of the apparent Coulomb-blockade dia- 
mond (Fig. [3| , a behavior which is actually observed in 
nanoelectromechanical systems. 



Planck equation ([6| . Numerical results for the scaled lin- 
ear conductance G = RI /V are shown in Fig.|3^-h, using 
the same parameters as in Fig. |3j:,f. We observe that the 
fluctuations reduce the size of the blockaded region and 
blur the edges of the conducting regions as the system 
can explore more conducting states in phase space. Nev- 
ertheless, the conclusions of the stability analysis clearly 
remain valid qualitatively. 



IV. LANGEVIN DYNAMICS 

To investigate the robustness of the stability analysis 
against fluctuations, we turn to the complete vibronic 
Langevin dynamics x -f ^{x)x = fcsix) + £,(t). The 
fluctuating force ^(t) is generated by fluctuations of the 
electronic occupation and the frictional force —^{x)x by 
the delayed response of the electrons to the vibronic 
dynamics. To compute 7(x) and £,{t), we solve Eq. 
([3| including the vibronic dynamics and the Langevin 
sources. Writing separate equations for average and 
fluctuations of the occupation by setting n — n + 5n, 
we see that the leading correction to h arises from the 
Poisson bracket, yielding n — uq — p Xdxno- At 
the same time, the fluctuations dn obey the correla- 
tor {6n{t)6n{t')) = r++r ^o(-^ ~ n-o)S{t — t'). Insert- 
ing these results into the expression for the current- 
induced force —gXn and employing reduced units, we 
find {£,{t)£,{t')) = D{x)5{t—t') with diffusion and damp- 
ing coefficients D{x) — 2ra;^no(l — rio)/v and j{x) — 
—rxd^no/v, respectively. Finally, we can pass fro m the 
Langevin to the equivalent Fokker-Planck equatiorP2l2ll 
for the probability V{x,p, r) that the nanobeam is at po- 
sition X and momentum p at time r, 

drV = -pd^V - f^tfdpV + jdpipV) + ^dlV. (6) 

Note that the diffusion and damping coefficients are non- 
vanishing only in the conducting region.l^^ 

The current / = ^ dxdpVst{x,p)I{x) is now obtained 
from the stationary solution drVst = of the Fokker- 



V. CONCLUSION 

We have presented a general approach to the interplay 
between continuous mechanical instabilities and current 
flow in nanoelectromechanical systems, and have applied 
our general framework to the Euler buckling instabil- 
ity. The current flow modifles the nature of the buck- 
ling instability from a continuous to a tricritical transi- 
tion. Likewise, the instability induces a novel tricritical 
current blockade at low bias. Our nonequilibrium Born- 
Oppenheimer approach generalizes not only to other con- 
tinuous mechanical instabilities, but also to other systems 
such as semiconductor quantum dots or single-molecule 
junctions with a discrete electronic spectrum, to other 
types of electron- vibron coupling,'2Sl and to further trans- 
port characteristics (e.g., current noise). 

Our proposed setup can be realized experimentally by 
clamping, e.g., a suspended carbon nanotube and apply- 
ing a force to atomic precision either using a break junc- 
tion or an atomic force microscope. Indeed, several recent 
experiments show that the electron-vibron coupling is 
surprisin gly str ong in suspended carbon nanotube quan- 
tum dotsPEGl 
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